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Abstract: We will give algorithms of computing bases of logarithmic cohomology 
groups for square- free polynomials in two variables. 

1 Introduction 

Let us denote by i? = C[x] — C[a:;i, . . . , a;„] the polynomial ring, by An = 
C{xi, . . . , Xn, di, . . . , dn) the Weyl algebra of order n over the complex numbers 
C and by (f^^, d) the complex of polynomial (or regular) differential forms (i.e. 
the complex of differential forms with polynomial coefficients) where d is the 
exterior derivative. 

The elements of An are called linear differential operators with polynomial 
coefficients. An element P{x, d) in An can be written as a finite sum P(a;, i9) = 
Eaaa(a;)5" where a = (ai,...,a„) € N", a^ix) € R and 9" = d^^ ■ ■ ■ d^" . 
Here di stands for the partial derivative 

For a non zero polynomial f ^ R we denote by i?^ the ring of rational 
functions 

and by (fi*, d) :— {Rj ®/j fi^, d) the complex of rational differential forms with 
coefficients in i?/ where d is the corresponding exterior derivative. 

Let us denote by Derc{R) the free i?-module of polynomial vector fields (or 
equivalently of C- linear derivations of R). Following K. Saito jT7] we will denote 
by DerR{— log/) the i?-module of logarithmic vector fields with respect to /, 
i.e. 

n 

Derni- log /) = {5 = ^ a^{x)d, G Derc{R) \ S{f) GR-f}. 

i=l 

Der{— log /) is canonically isomorphic to the _R-module Syzfi{di{f), . . . , dn{f), /) 
of syzygies among (9i(/), . . . , dn{f), /). This isomorphism associates the log- 
arithmic vector field S = J^i 0'i{x)di with the syzygy (ai(x), . . . , a„(a;), ~^y^). 
We will denote simply Der{— log/) if no confusion is possible. 

If / is a non zero constant, then Der{— log /) = Derc{R)- So we will assume 
from now that / is a non constant polynomial in R. 



It is clear that 

fDerciR) C Derni-logf) C Derc{R) 

and then Der{— log/) has rank n as i?"module. The i?-niodule Derni— log/) 
does not depend on the polynomial / but only on the hypersurface D = V(/) := 
{a e C" I /(a) = 0} C C". 

Assume / is reduced (i.e. / is square-free). According to K.Saito [T7] a 
rational differential p-form a> g fJ^ is said to be logarithmic with respect to / 
(or with respect to the hypersurface D = V{f) C C") if both fuj and fduj are 
regular (i.e. fuj e and fduj € ^^^^^)- We denote by ^^i\ogf) the i?-module 
of logarithmic differential p-forms with respect to /. K. Saito [TTl Corollary 
1.6] proved that Derfi{— log/) is a reflexive i?-module whose dual is r^^(log/). 
We denote by (f2*(log /), d) the complex 

r!0(log/) ^ r!^(log/) ^ . . . ^ f^"(iog/) — > 

which will be called the logarithmic de Rham complex and is also, for simple 
notation, denoted by SI* (log/) if no confusion arises. 

Algorithms of computing dimensions and bases of the de Rham cohomology 
groups H^{Vl') are given by T.Oaku and N.Takayama [T3], [IS] and U.Walther 
|19| . Here, / is any non-zero polynomial in n- variables. The purpose of this 
paper is to give algorithms of computing dimensions and bases of the logarithmic 
de Rham cohomology groups {VI' (log /) ) as C- vector spaces in the case of two 
variables. 

1.1 Logarithmic Comparison Theorem 

The rings R and i?/ have natural structures of left A„-module where di acts 
on a polynomial g and on a rational function -p^ as the partial derivative with 
respect to Xi. 

The de Rham complex of a left A„-module M, denote by DR{M), is by 
definition the complex of C-vector spaces (M (gjn Sl^, V*) where 

VP : M (g)R % M (g)R 

is defined, for p > 1, by yP{m ^ ui) — V"(to) A uj + m (E) dco and V°(m) = 
J2i diini) (8) dxi. Note that am (g) w = m aw for m G M, w € fi^ and a € R. 
The complexes fi* and DR{Rf) are naturally isomorphic. 

For any non zero f £ R, the inclusion if is a natural morphism of complexes 

tf:n'{\ogf)^n}. 

We say (see |3]) that / satisfies the (global) logarithmic Comparison Theorem 
if the morphism if is a quasi-isomorphism (i.e. if «/ induces an isomorphism 
HP{n'{logf)) ^ HP{n}) for any p). 

If 71 = 2, by [3l Cor. 2.7] and [2 Th. 1.3], if is a quasi-isomorphism if and 
only if / is a quasi-homogeneous polynomial. 
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1.2 The case n = 2. Bases for DerR{—\og f) 

If n = 2, any reflexive i?-inodule is projective and then, by Quillen-Suslin 
theorem, this i?-module is free. So, if n = 2, the i?-module Derii{— log f) is 
free of rank 2. In this case, we would like to compute a basis of Derji{logf) by 
taking the polynomial / = f{x,y) as input. By using the isomorphism 

Deri-logf) ~ SysR{di{f),d2if)J) 

and using Groebner basis computation, a system of generators of Derji[— log /) 
can be calculated. Then we can apply Quillen-Suslin algorithm (as presented 
for example in [5] and implemented in [B]) to compute such a basis. Known 
Quillen-Suslin algorithms use Groebner bases computation. Nevertheless, in 
some cases, for a big family of polynomials f{xi, X2) we will use an easier way 
to compute a basis of Der{— log/). 

First of all, we can assume / to be a reduced polynomial since Der[— log/) 
depends only on the affine plane curve D = V(/) = {(ai, 02) G | /(ai, a2) = 
0} C C2. 

Assume the plane curve D — V(/) is not smooth. The singular points of the 
plane curve D = V(/) (i.e. the affine algebraic set 

Szng{D) := V(/, /i, /a) = {a= (ai, 02) G | /(a) = /i(a) = /2(a) = 0}) 

-where /i = di{f), = c^2(/) ^ consists in a finite number of points (and it is 
not the empty set). 

We will consider the affine plane as a Zariski open subset of the projective 
plane P2(C), the affine point (ai, a-i) is mapped into the point with homogeneous 
coordinates (1 : ai : 02). Coordinates in P2(C) will be denoted by (xq : x\ : X2) 
and then the line at infinity is defined by xq — 0. 

Let us denote h = H{f),hi = H{fi) and ft,2 = H{f2) where H{—) denotes 
dehomogenization with respect to the variable xq. We will denote hy Z = 
Vp(/i, /ii, /i2) C P2(C) (resp. Z' ~ V(/i, ft-i, ft.2) C C'^ ) the projective algebraic 
set (resp. the affine algebraic set) defined by the polynomials h,hi,h2- The 
non-empty set Z (resp. Z') consists of a finite number of points in P2(C) (resp. 
a finite number of straight lines in C'^). Denote by 5 = C[xo, a;i, ^2] the ring 
of polynomials graded by the degree of the polynomials. If J = {h,hi,h2) 
denotes the ideal in S generated by h,hi,h2 then the quotient ring S/J has 
KruU dimension 1. Let us denote by 5*+ the irrelevant ideal in S, i.e. the ideal 
generated by xq,xi,X2. 

Proposition 1.1 The graded ring S/J is Cohen- Macaulay if and only if J is 
unmixed (i.e. is not an embedded prime associated with J). 

Proof: If S/J is Cohen-Macaulay then J is unmixed (see [5]). If J is un- 
mixed then S-i- is not an embedded prime of J and then the set of non zero- 
divisors of S/ J contains homogeneous elements of positive degree. That proves 
depth{S/J) > 1 but we also have depth(S'/ J) < dim(S'/ J) = 1. [] 
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If S/J is Cohen-Macaulay then the projective dimension of S/ J is 2 and J 
satisfies the Hilbert-Burch Theorem [5], i.e. there exists an exact sequence 



where (piigoi 91t 92) — goh + gi/ii + .92^2 and 02 is defined by a syzygy matrix 
of 4>i. In particular, since ker((/)i) = Syzs{h, hi, ^,2) is a graded free S'-module 
of rank 2 we can compute {s*^^-' = (sio, sn, S12), s'^^ = (s20, S21, S22)} a minimal 
system of generators and this system is in fact a basis of ker((/)i). By deho- 
mogenization (i.e. by setting xq = 1), we obtain a system '^[fo=iJ' 
generators of Syzji{f,fi,f2) — Der^X— log f) and since this i?-module is free 
of rank 2, this last system is in fact a basis. 

If S/ J is not Cohen-Macaulay we cannot apply, in general, the Hilbert-Burch 
theorem and the previous procedure fails to compute a basis of Dern^— log/). 

Example 1.2 (a) Consider the polynomial f — {x^ + xy^){x^ — j/^). With 
the notations as before (and writing xi = x,X2 = y, xq = t) we can use Macaulay 
2 to prove that the corresponding S/J is Cohen-Macaulay and to compute a min- 
imal system of generators of Syzs{h, /ii, /12) and then a basis of Derji{~ log/). 

Macaulay 2, version 0.9.2 

— Copyright 1993-2001, D. R. Grayson and M. E. Stillman 

— Singular-Factory 1.3b, copyright 1993-2001, G.-M. Greuel, et al. 

— Slngular-Libf ac 0.3.2, copyright 1996-2001, M. Messollen 

11 : R=ClQ[t,x,y] ; 

12 : f=(x"3+y'"4+x*y"3)*(x"2-y"2) ; 

13 : f l=dlf f (x, f ) ,f 2=dlf f (y , f ) ,h=homogenlze(f ,t) ,hl=homogenlze(f 1 , t) ,h2=homogenlze(f 2 , t) ; 

14 : Jf=ldeal(h,hl,h2) ; 

04 : Ideal of R 

15 : pdlm coker gens Jf 

05 = 2 

16 : Syzf=kernel matrlx({{hl ,h2,h}}) 

06 = Image {5} I x3+l/3x2y-4/3xy2 -tx2+4txy+3x2y+4xy2-y3 I 

{5} I 2/3x2y+l/3xy2-y3 tx2-txy+3ty2+2xy2+4y3 I 
{6} I -5x2-5/3xy+6y2 5tx-18ty-15xy-23y2 I 

3 

06 : R-module, submodule of R 

17 : mlngens Syzf 

07 = {5} I x3+l/3x2y-4/3xy2 -tx2+4txy+3x2y+4xy2-y3 I 

{5} I 2/3x2y+l/3xy2-y3 tx2-txy+3ty2+2xy2+4y3 I 
{6} I -5x2-5/3xy+6y2 5tx-18ty-15xy-23y2 I 
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3 



2 



o7 



Matrix R < R 




(b) Consider the polynomial g ~ {x^ + y** + xy^){x'^ + y^). With the notations 
as before (and writing xi — x,X2 = y, xq — t) we can use Macaulay 2 to prove 
that the corresponding S/ J is not Cohen-Macaulay and the minimal number of 
generators of Syzg{h, hi, /12) is 3. We can continue the last Macaulay 2 session: 

iS : g=(x"3+y"4+x*y"3)*(x"2+y"2) ; 

i9 : gl=dif f (x,g) ,g2=dif f (y ,g) ,h=homogenize(g,t) ,hl=homogenize(gl , t) ,h2=homogenize(g2 , t) ; 

110 : Jg=ideal(h,hl,h2) ; 

111 : pdlm coker gens Jf 

011 = 3 

112 : Syzg=kernel matrix({{hl,h2,h}}) 

012 = 

image 

{5} I tx2-5x3-4txy-20/3x2y-2xy2-5/3y3 x4+4/3x3y+x2y2+4/3xy3 tx3-tx2y+4x3y+4txy2+16/3x2y2+2xy3+4/3y4 I 

{5} I tx2+txy-10/3x2y-3ty2-5xy2-l/3y3 2/3x3y+x2y2+2/3xy3+y4 -txy2+8/3x2y2+3ty3+4xy3+2/3y4 I 

{6} I -5tx+25x2+18ty+100/3xy+ll/3y2 -5x3-20/3x2y-13/3xy2-6y3 -5tx2+5txy-20x2y-18ty2-80/3xy2-16/3y3 I 



012 : R-module, submodule of R 
il3 : mingens Syzg 

013 = 

{5} I tx2-5x3-4txy-20/3x2y-2xy2-5/3y3 x4+4/3x3y+x2y2+4/3xy3 tx3-tx2y+4x3y+4txy2+16/3x2y2+2xy3+4/3y4 I 

{5} I tx2+txy-10/3x2y-3ty2-5xy2-l/3y3 2/3x3y+x2y2+2/3xy3+y4 -txy2+8/3x2y2+3ty3+4xy3+2/3y4 I 

{6} I -5tx+25x2+18ty+100/3xy+ll/3y2 -5x3-20/3x2y-13/3xy2-6y3 -5tx2+5txy-20x2y-18ty2-80/3xy2-16/3y3 I 
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3 



ol3 



Matrix R < R 



We will revisit this example in Example 14.11 



2 Logarithmic 74„— modules 

Let us denote by Af'°s/ the quotient yl„~module M^°sf = 
Moreover, we denote by DerR{— log/) the set 



A„DerR(- log/) • 



Derni-logf) = {S+^\5e Dern{- log f)} 
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and by the quotient ^„-niodule 

M^ogf ^ ^" 

A^Derni-logf)' 

As quoted in subsection 11.21 for n = 2 the i?-niodule Der{— log f) (and 
hence il^{\ogf)) is free of rank 2. Moreover, by Tfl, 1.8] there exists a i?-basis 
{61,62} of Der(— log/) satisfying det(^) = / where 

6i = audi + aad2, i = 1, 2 

and A is the matrix (a^j). Then the dual basis of {61,62} is {(^1,072} with 

071 = — {a22dxi — a2idx2) uj2 = j{-ai2dxi + 011^x2). 

The i?-niodule ri^(log /) is free of rank 1 and ujiAuj2 is a basis of it. Moreover 
we have wi A ^2 = 

Proposition 2.1 Lei f € R — C[x,y] be a non zero reduced polynomial. There 
exists a natural quasi-isomorphism 

n' (log /) ^ nHoniA, (m'°s / ^ ^) 

where the last complex is the solution complex of ^1^°^^ with values in R. 

This Proposition is proven in • Ij in a more general setting using the notion 
of Vg-module. We will give here a direct proof to apply for our algorithm of 
computing logarithmic cohomology groups. 

Proof: F.J. Calderon [1] defines the so called logarithmic Spencer complex asso- 
ciated with AI^°^^ . In our situation, once a basis {^i, 62} is fixed in Der{~ log /), 
this complex is nothing but 

O^A^A^^A^O (1) 

where A stands for A2, the A-module morphism ci is defined by ei(Pi,P2) = 
Pi^i + P262 (for Pi e A) and £2 is defined by e2{Q) = Q(-^2 - bi,6i - 62) for 
Q € A and the polynomials hi being defined by the equality [^i,(52] = 6162 — 
6261 = bi6i + 62^2- In [I] it is proven that this complex is a A~free resolution 
of the module M^°sf_ We will use this resolution to find a complex of C-vector 
spaces representing the solution complex R7Jom^(M'°s , i?). Applying the 
functor HomA{—, R) to the logarithmic Spencer complex and using the natural 
isomorphism R ~ Ham a {A, R) , we obtain the complex 

O^R^R'^^R^O 

where e^g) = ((5i(g), ^2(3)) for g e R and €2{hi, /12) = ^1(^2) - 62{hi) - bihi - 
62 /i2 for hi e R. There is a natural morphism of complexes 
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r!0(iog/) = i? 4 n\\ogf) 4 n^iiogf) 

where r]o = id, r]i{hiuji + /i2'^2) = {hi,h2) and 772(3^^1 A ti;2) = for hi, h2, g G 
R and where {wi,tt'2} is the dual basis in il^(log/) of the basis {61,62} in 
Der{— log/). It is obvious that this morphism 77, of complexes of vector spaces 
is in fact an isomorphism of complexes. That proves the proposition. [| 

To each finitely generated left A„-module M we associate the complex of 
finitely generated right A„-modules RHomA^iM, An). To this one we associate 
the complex of finitely generated left A„-modules Homni^^, HHomA„ {M, An)) 
which is by definition the dual M* of the left A„-module M. 

If M is holonomic (i.e. if the dimension of the characteristic variety of 
M is n) then it can be shown that Ext\ (M, An) — for i n and then 
M* is the left holonomic A„-modulc HorriR{^V}^,Ext\^{M,An)) (see e.g. [TOl 
pag. 41]). Assume Ext\^{M,An) — ^ for some right ideal J C A„. Then 
HomR{il^, An/ J) is naturally isomorphic to the left A„-module jr where 
is the left ideal = {P^ \ P G J} and P^ is the formal adjoint of the operator 
P. 

If iVi , N2 are finitely generated left A„-modules there exists a natural iso- 
morphism of complexes 

RHoniA^ iNi,N2 )-^RHomA^ {RHoruA^ {N2,An), RHotjia^ (iVi , A„ ) ) 

and then a natural isomorphism 

RHorriA^iNi, N2)->RHomAjN; , N^). 

In particular, if N2 = R = C[a::i, . . . , a;„] then there exists a natural isomor- 
phism from RHoniA^iNi, R) (i.e. the solution complex of A'^i) to 

RHomAAR*,N^)- 

As the complex RHomA„{R, An) is naturally isomorphic to fl^ we can identify 
R and R* and then we have a natural isomorphism 

RHomA^{Ni,R) ^ RHomA^{R,Nl) ^ DR{N*). (2) 

Proposition 2.2 Let f € C[x,y] be a non zero reduced polynomial. Then there 
exists a natural isomorphism 

(A'&sf)* ~ M^°sf. 

Proof: This is one of the main results in 4]. We include here its proof for 
the sake of completeness. First of all, both A2-modules M'°s/ and M^°'^-f are 
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holonomic. That can be deduced from (TJ Cor. 4.2.2] since the set of principal 
symbols {a{Si), a{62)} is a regular sequence in the polynomial ring ^2] and 
then the KruU dimension of the quotient ring 

is 2. Then the characteristic variety of both yl2-modules and A'P"^^ has 

dimension 2 and the modules are holonomic. 

We will use the logarithmic Spencer complex associated with Af'°s f (see the 
complex ([1])) in order to compute Ext\{M, A) where A = A2 and M = M^°sf. 
Applying the functor HomA2{—, A2) to the complex ([1]) we get (by using the 
natural isomorphism HomA2{A2, A2) ~ A2) 

0~^A^A^^A^O 

where eT(P) = {5iP,S2P) and e^{Pi,P2) = {-S2 - + {5i - 62)^^2. Then 

we have 

Ext\{M, A) 



So, 



M* 



{-52 - 61, 5i - b2)A' 
A 



A{{-S2-hr,{6i-b2V)- 
Finally, {~62-bi)^ = 62 + ^ and (61-62)^ = -^i - (see H Cor. 3.1]). 



Theorem 2.3 For any non zero reduced polynomial f € C[a;,y], the complexes 
n'{\ogf) and DR{M^°^f) are naturally quasi-isomorphic. 

As a consequence of this theorem and by [13] , [13] and [TH] , the cohomology 
of the complex r2*(log /) can be computed starting with the given polynomial /, 

since a system of generators of the i?-modulc Derii{— log f) can be computed 
using the i?-syzygies of (5i(/), 92(/), /)■ _ _ 

Proof: Let us simply denote R = C[x,y], A = A2, M = A'f°^f , M = M^°^f. 
By Proposition 12.11 there exists a natural isomorphism 

17' (log/) ^ RHomA{M,R) 
and by equation ^ there exists a natural isomorphism 

RHomAiM,R) ^ DR{M*). 

By Proposition H^U we have DR{M*) ~ DR(M). 

We can give the explicit form of this quasi-isomorphism of complexes r* : 
^*{\ogf)^ DR{M). _ _ 

r° : i? ^ M is defined by T"(g) = gf where ( ) means the equivalent class 
in the ideal A2Derii{\ogf). 
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ri : f^Hlog/) M (g)R n\ is defined by 

i 

: Q'^{\ogf) M (^R nj^ is defined by T^{guJi A CJ2) g (81 /t^i A W2- Q 

3 Algorithm 

Let us summarize our algorithm of computing logarithmic cohomology groups 
in the two dimensional case. Most tensor products ® in the sequel are over A2. 
If we omit the subscript A2 for (8), it means that the tensor product is over A2. 

Algorithm 3.1 

Input: a non zero reduced polynomial f{x,y) 
Output: dimensions and bases of iJ'(f7*(log /)). 

1. Compute a free basis s = {sq, 81,82) and t = (toi^ii^2) of the syzygy 
module of f,fx,fy over the polynomial ring C[x,y\. This step can be 
performed by the following way. 

(a) Compute the minimal syzygy of h{f), h{fx), h{fy). Here, h{g) is 
the homogenization of 5. If the number of generators is 2, then the 
dehomogenizations of these generators are s and t. 

(b) If we fail on the first step, apply an algorithm for the Quillen-Suslin 
theorem to obtain s and t (call the procedure Quillen-Suslin). 

2. Define a left ideal in A2 by 

/ = ^2 • {-So + sidx + 82dy, -to + hdx + <29.y} . (3) 

Conipute the dimensions and bases of the de Rham cohomology groups 
for M = A2// with the algorithm in [13], [15]. In other words, replace the 
A2-module C[x,y, 1//] by A2/I of Q in the algorithm 1.2 in [H] , 

3. The bases of the previous step are given in A2/{dxA2 + dyA2)®M' where 
M' is (1, 1, —1, — l)-adaptive free resolution of M . Bases of de Rham co- 
homology groups in ri*®M ^q.i.s DR{M) ~q.i.s ^^*(log/) are determined 
by the transfer algorithm of U.Walther [13 Theorem 2.5 (Transfer The- 
orem)] and the correspondence t' given in our Theorem 12.31 Here, SI* is 
the Koszul resolution of the right yl2-module A2/{dxA2 -I- dyA2). 

In the first step, we should firstly try to find the minimal syzygy. Because, 
mostly it is faster than applying implementations and algorithms for the Quillen- 
Suslin theorem. 

The following example will illustrate how our algorithm works. 
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Example 3.2 We consider the case of / = xy{x — y). Two canonical generators 
of / = Z3erfl (log/) are 

^1=3 + xdx + ydy, £2 = -{2x - y) + (-.t^ + xy)dx 

The associated canonical logarithmic forms are 

uJi = jx{x - y)dy, L02 = j{~ydx + xdy) 

Let us proceed on the step 2. We apply the procedure of computing the de 
Rham cohomology groups [13], [H] for A2II. The maximal integral root of the 
6 function for I = A2- {-^i, -^2} with respect to the weight (1,1,— 1,-1) is 1. The 
dehomogenization of the (1, 1, —1, — l)-minimal filtered free resolution of A2II 
is 

A* : A2 [0] A2 [1] ® A2 [0] A2[l] (4) 

where 

a-^{c) = c(-^2,^i-l) fore 6^2 
a-\c, d) = (c, ) for (c, d) e A2 [1] ® A2 [0] 

Following [13,, procedure 1.8], we truncate the complex A2/ {dxA2+dyA2)^A2A* 
to the forms of (1, 1, — 1, — l)-degree at most 1 since the maximal integral root 
of the 6-function is 1 . The truncated complex is the following complex of finite 
dimensional vector spaces 

C^^{C + Cx + Cy)®C^^{C + Cx + Cy)^^0 (5) 

Here, 

d-\l) = {-£2,^1-1) mod d,A2 + dyA2 
= (0,0) 

d~^{a + bx + cy,d) = {a + bx + cy)£i + d£2 moddxA2 + dyA2 
= a 

Therefore, the cohomology groups H'^{A2/ {dxA2 + dyA2) A*) are 

H''{A2/{dxA2 + dyA2)®A') = Kera-^^C 

H\A2/{dxA2 + dyA2)<»A') ^ Kerd~^ /Imd-'^ ^ {Cx + Cy) ® C 

H\A2/idxA2 + dyA2)®A') = KeraVlma-i = Cx + Cy 

Finally, we perform the step 3. Put M — A2//. In order to give bases of the 
cohomology groups in M (^^ fi^, we apply the transfer theorem (algorithm) of 
Uh Walther [H]. 
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We consider the following double complex (c.f., 2.4 of T^). 



n{2) (g) A2 n{2) (g) {A2 e A2) 



•A2 



A2 



(g) ^2 



A2 ® {A2 ® A2) 



A2 

A2 



A2 ® {A2 ® A2) 



■f7(2) 



^2 ® A2 



A2 



A. 



A2 «) A2 ■ 



■ ^2 «) M 



A2 



A2 



M 



■ y42 (g) M 



Here we denote A2l(dxA2 + dyA2) by ri(2), which is isomorphic to fi^ as the 
right A2-module. The vertical complex is constructed by the Koszul resolution of 
r2(2) as the right module denoted by il' . The horizontal complex is constructed 
by A'. Note that we have the following maps in the complex: 



2((a,6)(g)c) = 



{—dyQ + dxb) ® c 
{dxCL, dyo) (g) c 



-1,-1 



((a,6)(g(c,d)) = 



-O.-l 



(a(g (c, d)) 



{dya - dxb) (g (c, d) 
{dxCL, dyO) (g (c, d) 



e^'"((a,6)®c) 



{—dya + g) c 
{dxCL, dyO) g) c 



^(ag)c) 
"\a(g (c, d)) 



a(g c(-f2,4 - 1) 
a(g (c^i +d^2) 



a^'"^((a,6) (gc) 
a^'"^((a,5) (g (c, d)) 



(a,5)(gc(-^2,^i-l) 
(a,6) (g (c^i + d^2) 



t°'-2(a(gc) 
^(a(g (c, d)) 



a(g c(-4,4 - 1) 
a (g (c^i + di2) 



The last vertical complex is quasi isomorphic to DR[M). Let us compute trans- 
fers. Two cohomology classes x and y in Kera° C Vt{2) (g) A2 are lifted to 1 ® a; 
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and 1 (g) y in A2 (8 ^2 respectively, and we push them to M. It follows from 
the definition of r^, xuji A UJ2 and yuji A UJ2 is the basis of i?^(f2(log /)*). 

Let us compute transfers of bases of iJ^(Sl(2) (g) A*). The cohomology class 
1 (g) (x, 0) in Ker are lifted to 1 (g (x, 0) in A2 g) (A2 © A2). We have a2'"^(l g) 
(a;,0)) = 1 (g x-^i. Solving —dyU + 9^:6 = cc^i in A2, we obtain the preimage 

MA ^ 

by 6"'^''^; we have e^'^{{—xy,x'^) g) 1) = cc^i. Push this clement to © g) M, 

we obtain (xydx — x^dy) g) 1. Let us compute the preimage by t^. Solving 
ci/wi + C2/W2 = xydx — x^dy, we obtain ci = 0, C2 = — x. Therefore, 1 g) (x, 0) 
stands for — xct'2. Analogously, 1 g) (y,0) is transfered to —y'^dx + xydy and 
stands for — ya;2 and 1 g) (0, 1) is transfered to x(y — x)dy and stands for wi. In 
summary, 

H\ni\ogf)') = C(-xV2 + C(-yV2 + Cuji. 
Finally, we compute transfers of bases of H^{n{2)^A*). Since a^'~^(lg)l) = 
1 gi {—£2,(1 — 1), we firstly need to compute the preimage of this element by 
e^'~^. Since the projection of this element to Q{2) g) {A2 © A2) is zero, we have 
—£2 = dxx{x — y) and €1 — 1 = dxX + dyy. We decompose 1 g) (—£2, — 1) as 

l(g (-^2,0) + Ig) (0,£i - 1) = -£2<»(l,0) + (4-l)<»(0,l) 

= dxx{x - y) ®{l,Q) + {dxX + dyy) ®{Q,l) 

Since e^'~^ is linear, this sum is equal to e^'^^{c) where ci = (0, — x(x — y)) g) 
(1,0) + (y, -x) g) (0,1). Since a^-'^{ci) = iyi2,-x{x ~ y)ii - XI2) ®l = 
{dxxy{y — x),dyxy{y — x)) g) 1, the preimage of a^''^{ci) by e°'° is equal to 
xy(y — x)(E)l G A2 (g M. Therefore, the preimage of r° is equal to —1 and hence 
7J°(0(log /)*) — C(— 1). Although we have done this computation by hand, 
computation of transfers can be done by Grobner basis computation. See [19| 
and the source code for deRhamAll of the Macaulay 2 package for D-moduIes 

Before presenting implementations and larger examples, we explain a bit 
about a procedure to find a preimage of r* in general. The transfer algorithm 
gives an element in gj^a M where ft' is the Koszul resolution of 11(2) ~ $1^ 
as the right A2-module. This element can be identified with a differential form 
with coefficients in M and we need to find the preimage of it by which lies in 
fi'(log/). This can be performed by the method of undetermined coefficients. 

Consider the case of r^. Take an element cilui + C2LO2 in f2^(log/) where 
Ci £ R. We have seen in Theorem 12.31 that 



^2 _ 

T^idLOi + C2UJ2) = fL0i'SlA2 Cl + fuJ2<»A2 C2 £ I © \ <Sl A2 M (6) 



A 



2 , 

Here, we identify ^ q ^ g) A2 '^^i with mi g)fl dx and (^^^ '^A2 "^2 with m2 ®iidy, 
rrii € M (comparison theorem) and when Ui = aidx + bidy, we denote ( ^* j 
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uJi. As the output of the transfer algorithm, we are given an element midx + 
m2dy, rrii G M . We regard rrii as an element in A2 in the sequel. We rewrite 
fijji as /wi = Adx + Bdy and fuj2 — Cdx + Ddy. Assume / is generated by ii 
and £2- Then, the definition of 1^ induces the following identity in A2 by 
taking coefficients of dx and dy 

2 

Aci + Cc2 = rrii+Y^ ^i^J + C^) 
2 

Bci +Dc2 = W2 + ^ + (8) 
j=i 

where Ci € R, d^ , e £ A2 are unknown. Fix a degree bound m for these elements 
and determine these elements by the method of unknown coefficients. The iden- 
tities d?]) and ^ induce a system of linear equations over C for the coefficients. 
Increasing the degree bound and solving the system, we will be able to obtain 
ci and C2 in finite steps by virtue of Theorem 12.31 

Consider the case of r^. Since our basis in H^{Q.' (g) M) is given in terms of 
X and y and fuJi A UJ2 ^ dx A dy, we need no computation to find the preimage 
byr2. 

Let us consider the case of r°. Let m be an output of the transfer algorithm. 
It lies in A2 in general. Finding the preimage g of r*' can be done by solving 
gf = m + '^j^j where g (z R and dj e A2. 

4 Implementation and Examples 

The second and third steps of Algorithm 13.11 can be performed with the help of 
the D-module package on Macaulay2; use the commands DintegrationAll to 
obtain the dimension of the cohomology groups, DintegrationClasses to ob- 
tain the bases of cohomology groups, and a modification of DeRhamAll to obtain 
the bases of cohomology groups in fi* (g) M. Unfortunately, this implementation 
has not installed an efficient algorithm of computing 6-function by Noro [11] to 
get the truncated complex in [13], [15]. Then, only relatively small examples 
are feasible. The Example 14.11 is computed by our Macaulay2 program. The 
Example 1121 is computed by our implementation on kan/kO and Risa/Asir with 
an implementation of lllj (the transfer algorithm has not been implemented 
yet for kan/kO). This implementation also uses the minimal filtered resolution 
to reduce the size of complex of A2-modules [16]. The program is contained 
in the OpenXM package with the name logc2.k (http://www.openxm.org). 
Our implementation does not contain that for the Quillen-Suslin theorem. We 
utilize the implementation by A.Fabianska on Maple when the step l-(a) fails. 
We also note that computation of the preimage of may become a bottleneck 
of computation. 
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Example 4.1 (Continued from Example 11.21 (b).) We will determine bases of 
i7*(ri' (log/)) where / = (s^ + J/'* + xy^){x'^ + y"^). We firstly use Fabianska's 
program for the Quillen-Suslin theorem to find the 2 free generators of the 
syzygies of /, fx, fy The two rows of the following matrix S are the generators 

q ^ (Sii {-23/Gy + l/2)x^ + (y^ + - 2y)x - 5/6y3 (l/3y + l/2)x^ + (-3j/2 + l/2y)x + y^ + ^jZy^ - 3/2?;2 
V'S'21 -46/75x3 + {4/25y2 - 2/25y)x2 - S/lSy^x 4/75x3 - 12/25yx2 + (4/25i/3 - 2/751/2)0; - 2/51/3 

where = (115/6j/-5/2)j;-62;3-43/6y2 + 9?;, ^si = 46/15x2 + (-24/25y2 + 

22/75y)a; + 12/5y2. Put A = f f f . Then, det(A) = i/. We put 

wi = j{a22dx — a2idy) and W2 = j{—ai2dx + andy). {s/ioji agrees with the Wj 
in Theorem O) 

We apply the integration algorithm ajnd the transfer algorithm for M . We 
obtain the following result. (1) H^{DR{M)) is spanned by 1 ® / and then we 
have i?°(n*(log/)) ~ C- 1. (2) H^{DR{M)) is spanned l(g)a where a runs over 

3 3 2 3 4 

o9 = {{1}, {-x}, {y }, {-x*y }, {x*y }, {x y}, {y » 

(We have pasted the output of our Macaulay 2 program trans .m2.) Then, we 
have 

H^Q'ilog /)) ~ (C • 1 + C • {~x) + --- + C-y^)ujiAL02 

(3) H^(DR(M)) is spanned by 3 differential forms midx + m2dy where mi, TO2 
are elements in A2, of which explicit expressions are a little lengthy. We solve 
the identities ([7]) and ([8]) to find ci and C2- In other words, we need to compute 
preimages of midx + m2dy by r^. As we explained, this can be done by the 
method of undetermined coefficients degree by degree. We can find solutions 
when the degree of Ci,dl,e with respect to a;,y is 6 and that with respect to 
dx, dy is 0. Here is a basis of 3-dimensional vector space iJ^(ri' (log/)) obtained 
by this method. 

—yxuji — 4/25x2^^2 

((215/28y- 1101/280)a;- 367/56y2)wi + (43/35a;2 - 367/3502/x)w2 

((y - ll/30)a; - 28/9y3 - IS/Gy^ + 14/3?/)wi + (4/25a;2 + (-112/225y2 + 2/by)x + 56/45?/2)cj2 
All programs and session logs to find this answer is obtainable from 

|http : //www ■ math . kobe-u ■ ac . j p/OpenXM/Math /LogCohomology/2007 - ll/log-2007-11-22 . txt 

The logarithmic comparison theorem does not hold for this example. In fact, 
the dimensions of the de Rham cohomology groups iJ*(ri'), {i — 2,1,0) are 
5, 3, 1 respectively. 

Example 4.2 We apply a part of our algorithm to compute the dimensions of 
the cohomology groups H^{Q.'{\ogf )) for f = x^ + y'^ + xy''^^ . Here is a table 
of p, q and the dimensions of , , 7?" and timing data. 
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p 


q 


Dimensions 


Timing in seconds 


10 


11 


(8,1,1) 


3.5 


10 


12 


(9,1,1) 


4.6 


10 


13 


(10,1,1) 


6.9 


10 


14 


(11,1,1) 


9.4 


10 


20 


(17,1,1) 


55.0 


10 


21 


(18,1,1) 


86.8 



The program is executed on a machine with 2G RAM and Pentium III (IG Hz). 

The homogcnization of /, fx, fy generates an ideal that is Cohcn-Macaulay. 
These examples do not need to call the subprocedure Quillen-SusUn. However, 
the logarithmic comparison theorem does not hold for these examples. Compu- 
tation of de Rham cohomology groups is not feasible by our implementation. 

5 A Yet Another Algorithm 

In the previous section, we have presented a general algorithm of computing a 
basis of the logarithmic cohomology groups for plane curves. However, this al- 
gorithm relies on algorithms for the Quillen-Suslin theorem and they are some- 
times slow. We will present a yet another algorithm, which is free from the 
Quillen-Suslin theorem, but it works only for computing a basis of the middle 
dimensional cohomology group i/^(il*(log/)) under some conditions on /. This 
section can be read independently from other sections. For reader's convenience, 
we will also redefine some notations. 

Before stating the main algorithm, we start with an introductory example, 
which explains the idea of our algorithm. 

Put K = C and L = (1 - x)xd + 2x{= 0^ - x(6'^ - 2)). We consider the 
problem of determining a basis of the ii'- vector space K[x\/ L- K[x]. Since L is a 
ii'-linear map and K[x?\ is an infinite dimensional if-vector space, the quotient 
has the structure of a if- vector space. However, note that L ■ K[x] is not an 
ideal and we cannot use Grobner basis to get a basis. 

Let us act L on monomials; L ■ x^ = kx^ — (fc — 2)x'^+^. For small k, they 
are L- 1 = 2a;, L-x = x + x"^, L-x"^ = 2x'^. Then x''+^ ~ k^^'' modulo L-K[x]. 
In particular, if A; > 3, then the monomial can be reduced to a lower order 
monomial modulo L ■ K[x]. Hence, the set of monomials l,x,x'^,x'^ generates 
K[x]/L ■ K[x]. More precisely, we can prove that it is isomorphic to Fs/L- F2. 
Where Fk is the set of polynomials of which degree is less than or equal to k. 
The monomials l,x,x^,x^ are not independent modulo L ■ K[x] and satisfies 
the relation above. Finally, we conclude that i('[x]/LJr[a;] 2:^ K ■ 1 + K ■ x^. 

Note that 3 is the magic number, which is characterized as follows. Put 
L* = —{1 — x)xd — 1 + 4x. 111(^1 _i^{L*)nK[—dx] is generated by 6(— 9a:;) where 
6(s) = s — 3. The polynomial b{s) is called the indicial polynomial (6-function) 
for integration. The magic number 3 is the root of b{s) = 0. We will call 
the method to bound a degree by a root of a 6-function b-function criterion. 
T.Oaku firstly introduced the 6-function criterion to compute restrictions and 
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integrations of D-modules [12]. The topic of computing K[x\/L ■ K[x\ by the 
6-function was also discussed in more detail in an expository book "D-modules 
and Computational Mathematics" (in Japanese) by T.Oaku. 
Let / be a polynomial in two variables. Put 

fij = fc-form with coefficients in K[x^y, 1//] 

As we have explained in the introduction, the k form lo G $7^ is called log- 
arithmic fc-form iff both of Joj and df A lo have polynomial coefficients. The 
space of logarithmic fc- forms is denoted by il.''{\ogf). The question we ad- 
dress in this section is the computation of ^jyf|^9)' easy to see that 

ri^(log/) = ^ — . Let us determine all the logarithmic 1-forms. Let 

(p, g, r) a triple of polynomials such that 

fyP - fxq + fr = (syzygy equation). (9) 

Note that {0,f,fx), (/, 0,-/^), {fx,fy,0) are trivial solutions of the syzygy 
equation. For a solution {p, q, r) of the syzygy equation, lo = p'^^+i'^v belongs 
to f2^(log/). Conversely, any logarithmic 1-form can be expressed in this way. 
In fact, the condition that df A lo has a polynomial coefficient is equivalent to 
that fyP — fxQ is a. multiple of /. 

Put LO = Pd£+My_ Let eix,y) be any polynomial. Then, d{eLo) ^ {Le)^^^ 
where 

L = qdj: - pdy + qx-Py + j'^'^'^ 

We denote the Weyl algebra A2 by D for simplicity in the sequel. Sup- 
pose that Li, [i — l,...,m) stand for a set of generators of the solution 
space of the syzygy equation, which is a if[a;, ?/]-module. Then dVL^ {log f) = 
^ LiK[x, y]dxAdy / f . Therefore, the computation of is nothing but the com- 
putation ofK[x,y]/ 1 L,9K[x, y]. Put /*=£>• {L^, ... , L*^}, which is a left 
D ideal. We denote by the iC-subvector space of D of which (1, 1, —1, — 1)- 
order is less than or equal to fc 18, p. 14, p. 203] 



Algorithm 5.1 H"^ {Q- (log f)) . 

Step 1. Find generators of the syzygy equation and obtain explicit expressions 
of L^. 

Step 2. Compute (1, 1, —1, — 1)-Gr6bner basis (standard basis) of /. We denote 
the elements of the Grobner basis by U (renaming) . 

Step 3. Find the monic generator h{—dxX — dyy) of in(i ^ _i (/) n K[—dxX — 
dyy\. 

Step 4. Let fcg be the maximal non-negative root of 6(s) = 0. Then, return 
/C- vector space basis {ci} of 



Fkol X! ' -^feo-ordd. 1,-1,-1) (L, 
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{cidx A dy/f} is a basis of H"^. 

The steps 2, 3, 4 can also be done by computing D/{I* + d^D + dyD) (O-th 
integral module) where /* is the formal adjoint of /. (As to details for the steps 
2, 3, 4, see [H].) 

Note: Although our discussion is independent from the discussions of the pre- 
vious sections, the left ideal generated by L* is nothing but Derji{— log/) and 
hence this algorithm and the Algorithm 13.11 are analogous for computing a ba- 
sis of H^{Q'{f)). We also note that finding bases for W{n' (g)^^ M) can be 
performed by applying the integration algorithm and the transfer algorithm for 
D/I*. The Algorithm 13 . 1 1 relies on algorithms for the Quillen-Suslin theorem to 
find bases for ff*(f7(log /)'), i = 1,0. 

Theorem 5.2 // dim /,) < 0, dimy(/,/,) < 1, dimV{fJy) < 1, 

then the Alaorithm \5. 1\ is correct. 

Proof: Let / be the left ideal in D generated by Li, . . . , L,„. We may assume 
that / contains fdx, fdy and fydx — fxdy. Therefore, the characteristic variety of 
/ is contained in V{f{x, y)^, /(x, y)r], fy(x, y)£, — fx(x, y)r]), of which dimension is 
less than or equal to 2 from the assumption. In fact, assume (a, b) £ V{f, fx, fy). 
Then, ^ and 77 are free and then the dimension of the characteristic variety is 
less than or equal to 2. Assume (a, b) e V{f, fx) \ V{f, fx, fy)- Then, we have 
f{a,b) = 0,fx{a,b) = and fy{a,b) ^ 0. Then, r\ is free and C = a-nd then 
the dimension of the characteristic variety is less than or equal to 2. The rest 
cases can be shown analogously. Therefore, Z?// is a holonomic _D- module and 
hence a non-trivial b exists ([T51 Chapter 5, Theorem 5.1.2]). The rest of the 
correctness proof is analogous with that of the 0-th integration algorithm of 
IJ-modules [12], [H Chapter 5; Theorems 5.2.6 and 5.5.1]. [] 

Note: The algorithm works to get _ff"(r2*(log /)) in the n-variable case if 
dimy(/, fxi T ■ ■ 1 fxi^) < n — m for all m — 1, . . . , n and all combinations 
The algorithm and the correctness proof are analogous. In fact, 
since f^i and {—^Y fxj£,i — {~^)'' fxi£.j , (1 < * 7^ j < "-) are in the characteristic 
ideal for / and then the dimension of the characteristic variety is less than or 
equal to n by utilizing the condition. 

Example 5.3 FoTf = {x^ + y'^ + xy^){x^+y^),weheivedimH^{n'(\ogf)) = 7 
with our yet another algorithm 1 5. II The execution time is 1.9s. We need to call 
the procedure Quillen-Suslin if we use the first algorithm. 

We close this paper with a final note and the acknowledgement of this paper. 
We think that logarithmic differential forms give nice simple bases for some of 
hypergeometric integrals as pairings of twisted cycles and cocycles when the 
logarithmic comparison theorem holds for twisted de Rham complex. We hope 
that our result have applications to study hypergeometric integrals. The authors 
are grateful to A.Fabianska for helping us to compute free bases of syzygies by 
using her implementation for Quillen-Suslin's theorem. 
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